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Abstract. The wasp spider Argiope bruennichi (Scopoli, 1772) is of Mcditerranean-Pontian origin, but for decades it has 
been expanding northwards, including into the territory of Poland. Based on well-documented expansion records, we can 
distinguish “old” (south-eastern and south-western) and “new” populations (north-eastern), respectively, from the 1930s to 
the 2000s. In Poland, some populations of A. bruennichi were expected to be more genetically isolated from others, due to 
distance effects or differential times of arrival. We evaluated whether the oldest populations were in a state of Hardy- 
Weinberg equilibrium ( HWE ), and whether recently founded populations were in an expansion phase. Specimens of A. 
bruennichi (n = 184) were collected at six localities in Poland and single sampling sites in Italy and Japan. Nine 
microsatellite loci were amplified although only five were useful in the final analyses. Based on the genotypes obtained, we 
estimated basie measures of genetic diversity and tested for deviation from HIVE. The results showed a low level of 
polymorphism amongst the investigated markers, and accordingly, we found a low genetic diversity in populations. Only 
populations from Italy and Japan, and one population from Poland, were in HWE. The level of genetic differentiation 
among sampling sites from Poland was also very low. The high dispersal ability of the wasp spider appears to have 
facilitated high gene flow among populations. The peripheral and recently settled populations were characterized by the 
highest heterozygosity and the lowest inbreeding coefficient (F IS ) The remaining Polish populations arc therefore still in 
the expansion phase, as indicated by deviations from HWE. 
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The current distributions of species in Central Europe 
reflect the continent’s glacial and post-glacial history, as well 
as the impacts of human activity on the environment. Habitat 
changes, such as a reversion to early successional stages and an 
increase in fallow areas, favor the migration and development 
of terrestrial arthropod populations. Habitat changes may 
inhibit range expansion, especially if the distances between 
patches are too large (Grez et al. 2004). However, in the case 
of organisms with high dispersal abilities, movement is still 
possible over long distances (Nee & May 1992; Kareiva & 
Wennergren 1995). 

Spiderlings or even some adult spiders can spread through 
the air over large distances (sometimes even several thousand 
kilometers) by ballooning on lengths of silk (Gressitt 1965). 
This allows an unstable or newly accessible environment to be 
inhabited (Meijer 1977). Species with a high potential for 
dispersal generally display a homogeneous population struc¬ 
ture, even over considerable distances (Slatkin 1993). In the 
case of orb-weaving spiders, whose dispersal capabilities are 
usually very high, a similarly high level of gene flow can often 
be observed among populations (Ramirez & Fandino 1996; 
Ramirez & Haakonsen 1999; Lee et al. 2004; Jung et al. 2006). 
Many studies have investigated genetic differentiation among 
populations of spiders (e.g., Ramirez & Fandino 1996; 
Ramirez & Haakonsen 1999; Lee et al. 2004; Jung et al. 
2006; Croueher et al. 2011; Krehemwinkel et al. 2016b), but the 
balance between genetic drift and migration at the margins of 
a species’ range is still poorly understood. Krehenwinkel & 
Tautz (2013) and Krehenwinkel et al. (2015) examined the 
genetic variability in native and invasive populations of the 
“wasp spider” Argiope bruennichi (Scopoli, 1772), revealing 


marked distinctions between them. They also pointed to a 
significant role for genetic admixture, as reflected in a high 
genetic diversity in populations at the edges of the range. 

Argiope bruennichi is of Mediterranean-Pontian origin, and 
its current range extends over the Palearctic region (Jager 
2012). During the past decades, populations of the wasp spider 
have been expanding northwards. In Poland, the first localities 
appeared in the west (1930s) and in the south-east (1960s), but 
this species was still rare until the 1990s (Urbanski 1948; 
Barabasz & Gorz 1998). In the late 1990s, the number of 
localities at whieh the wasp spider increased considerably, and 
the species was soon present across the whole territory of 
Poland, except in the north-east where the climate is 
particularly harsh. Low temperature seems to matter most 
where the expansion of thermophilous species is concerned, 
and is considered as a limiting factor (Guttmann 1979; 
Kumschick et al. 2011) that probably explains the later arrival 
of A . bruennichi in the northern part of Europe (Jonsson 2004; 
Terhivuo et al. 2011). The rapid expansion in range may be the 
result of climatic warming, or else the local adaptation of the 
species to lower temperatures (Kumschick et al. 2011). Some 
studies have shown that an invasive population can exhibit 
increased cold tolerance, e.g., through the generation of an 
admixture of European and Asian lineages, the latter probably 
being pre-adapted for lower temperatures (Krehenwinkel & 
Tautz 2013; Krehenwinkel et al. 2015, 2016a). It is nevertheless 
likely that, in such newly-eoionized areas as the mountains of 
southern Poland or the north-east of Poland, populations of 
A. bruennichi will be vulnerable to severe temperature 
fluctuations to the point at which effective long-term 
colonization may be hindered. On the other hand, stable 
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Figure 1.—Map showing the sampling locations. 


populations can be expected to persist in the areas that were 
occupied first. 

It is now recognized that genetic diversity decreases with 
increasing distance from the source population, and this is also 
true of invasive species (Ciosi et al. 2011; Kim & Sappington 
2013; Kononov et al. 2016). Populations at the limits of a 
species’ range differ from those at the source in manifesting 
reduced levels of polymorphisms that reflect the overall loss of 
genetic variation. Likewise, at a finer geographic scale, the 
further from the center of the habitat, the lower the 
heterozygosity (Peter & Slatkin 2013). Hardy-Weinberg 
disequilibrium has been observed in distinct locations for 
certain invasive species, across several generations (Zima et al. 
2016). Moreover, edge populations can show low genetic 
variation when an investigated generation derives from just a 
few founder migrants (Bunn & Hughes 1997). We therefore 
hypothesized that due to the small sizes and potentially 
ephemeral nature of founding populations, the genetic 
diversity of A. bruennichi would be lowest in north-eastern 
Poland, where the climate is harsh and the species only 
appeared recently. We in turn expect that the oldest 
populations - from the west and south-east of Poland - are 


in a state of Hardy-Weinberg equilibrium and present higher 
genetie diversity, while populations from the north-east of 
Poland remain in the expansion phase and are thus more likely 
to deviate from Hardy-Weinberg equilibrium. 

METHODS 

Study sites.—Specimens of A. bruennichi were collected 
from six localities in Poland (Fig. 1), taking into account the 
time of colonization and the presence of dispersal barriers 
(e.g., forests, rivers, mountains and human structures). These 
localities were: SUW (harsh elimate conditions, youngest 
population, arrival ca. 2005, at 54°18'50"N, 22°27T8 // E); TOR 
(mid-forest meadow habitat, oldest population, arrival ea. 
1930s, at 52°14T8 // N, 15°06T6"E); WAR (eity habitat, at 
52°11 '30"N, 21°04'07"E); WOL (island habitat, at 
53°50'56"N, 14°19'43"E); WYS (open area habitat without 
barriers, at 50°56'54"N, 22°43'36"E); ZAW (mountain valley 
habitat, at 49°40 / 20 // N, 19°37'43"E). Additional samples came 
from “old” populations in Italy (1 locality, at 40°49 / N, 
16°25'E) and Japan (1 locality, at 34°15'N, 132°34 / E). We 
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Table 1.—Characterization of mierosatellite polymorphisms in eight populations (//= 184) of A. bruenniehi. n — sample size within population 
with successful amplification; A — number of alleles; R A — allelic richness; R ? — private allelic richness; E A — effective number of alleles; H 0 — 
heterozygosity observed; H E — heterozygosity expected; HIVE — P -values for HWE exact test for heterozygote deficiency/excess (those 
statistically significant arc in bold); F IS — fixation index;* — statistically significant values of F 1S and P -values for HWE after Bonferroni 
correction (800 randomizations, adjusted P -value at a = 0.05 was 0.0013). 


Locus 

A 

Ra 

Rp 

Ea 

Ho 

He 

HWE 

b'is 

MA3 

10 

7.31 

0.11 

2.84 

WAR (n = 25) 
0.480 

0.648 

<0.05 

0.278 

MA5 

7 

5.64 

0.02 

2.80 

0.680 

0.642 

0.616 

-0.038 

MA7 

2 

2.00 

0.00 

1.77 

0.400 

0.435 

0.663 

0.101 

MA33 

6 

5.34 

0.71 

3.58 

0.320 

0.721 

<0.001* 

0.570* 

MA35 

3 

2.81 

0.03 

1.44 

0.360 

0.306 

1.000 

-0.155 

MA3 

8 

7.00 

0.13 

4.48 

TOR (n = 25) 

0.320 

0.777 

<0.001* 

0.601* 

MA5 

5 

4.78 

0.08 

3.25 

0.720 

0.692 

0.751 

-0.020 

MA7 

3 

2.56 

0.26 

1.71 

0.320 

0.414 

0.484 

0.246 

MA33 

4 

3.97 

0.00 

2.69 

0.360 

0.628 

<0.01 

0.443 

MA35 

2 

2.00 

0.00 

1.32 

0.200 

0.241 

0.386 

0.189 

MA3 

6 

5.37 

0.00 

2.28 

WOL (i n = 25) 
0.423 

0.561 

<0.05 

0.265 

MA5 

6 

5.36 

0.00 

2.73 

0.654 

0.634 

0.757 

-0.012 

MA7 

3 

2.54 

0.24 

1.36 

0.308 

0.266 

1.000 

-0.140 

MA33 

3 

3.00 

0.00 

2.75 

0.154 

0.636 

<0.001* 

0.766* 

MA35 

2 

2.00 

0.00 

1.30 

0.269 

0.233 

1.000 

-0.136 

MA3 

10 

7.76 

0.57 

4.58 

WYS (n = 25) 

0.640 

0.782 

0.126 

0.201 

MA5 

8 

6.87 

0.07 

4.11 

0.600 

0.757 

<0.01 

0.227 

MA7 

3 

2.56 

0.56 

1.76 

0.520 

0.431 

0.730 

-0.186 

MA33 

5 

4.78 

0.15 

2.98 

0.240 

0.665 

<0.001 

0.651* 

MA35 

4 

3.12 

0.81 

1.56 

0.440 

0.358 

0.710 

-0.208 

MA3 

11 

8.66 

1.04 

4.39 

SUW (#i = 25) 
0.640 

0.772 

<0.05 

0.191 

MA5 

8 

6.88 

0.59 

4.22 

0.640 

0.763 

0.176 

0.181 

MA7 

2 

2.00 

0.00 

1.85 

0.640 

0.461 

0.091 

-0.371 

MA33 

4 

3.37 

0.00 

2.09 

0.520 

0.522 

0.525 

0.025 

MA35 

4 

3.37 

0.28 

2.09 

0.560 

0.522 

0.899 

-0.052 

MA3 

9 

7.52 

0.42 

4.03 

ZAW (n = 25) 

0.800 

0.752 

0.978 

-0.043 

MAS 

6 

5.10 

0.05 

2.58 

0.600 

0.613 

0.230 

0.041 

MA7 

2 

2.00 

0.00 

1.27 

0.160 

0.211 

0.287 

0.262 

MA33 

3 

3.00 

0.00 

2.68 

0.200 

0.626 

<0.001* 

0.692* 

MA35 

3 

2.81 

0.03 

1.56 

0.200 

0.358 

<0.05 

0.457 

MA3 

8 

6.71 

0.00 

2.52 

1TA (n = 20) 

0.700 

0.604 

0.341 

-0.134 

MAS 

4 

3.70 

0.00 

2.69 

0.550 

0.629 

0.668 

0.150 

MA7 

3 

2.62 

0.70 

1.16 

0.050 

0.141 

<0.05 

0.661* 

MA33 

2 

2.00 

0.00 

2.00 

0.550 

0.499 

1.000 

-0.077 

MA35 

2 

2.00 

0.00 

1.92 

0.300 

0.480 

0.153 

0.397 

MA3 

4 

4.00 

0.15 

2.97 

JAP (ii = 14) 

0.357 

0.663 

<0.05 

0.490 

MA5 

4 

4.00 

0.09 

2.35 

0.357 

0.574 

0.066 

0.409 

MA7 

3 

3.00 

1.00 

1.81 

0.357 

0.446 

0.316 

0.235 

MA33 

4 

4.00 

1.00 

1.46 

0.357 

0.314 

1.000 

-0.102 

MA35 

6 

6.00 

4.00 

3.73 

0.643 

0.732 

0.125 

0.158 
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Table 2.—Mean allelic diversity and heterozygosity indices for five microsatellite loci in eight populations (n = 184) of A. bruennichi. A — 
mean number of alleles; R — allelic richness; P A — private alleles; Rp — private allelic richness; E A — effective number of alleles; H Q — 
heterozygosity observed; // H — heterozygosity expected; HWE — P-valucs for HWE exact test for hetcrozygolc deficiency/excess; F iS — fixation 
index; * — statistically significant values of F iS and P-values for HIVE after Bonferroni correction (800 randomizations, adjusted P-value at a = 
0.05 was 0.0013). 


Population 

A 

R 

Pa 

R P 

Ea 

Ho 

He 

HWE 

^1S 

WAR 

5.60 

4.62 

0.20 

0.17 

2.49 

0.448 

0.551 

<0.001 

0.206 

TOR 

4.40 

4.06 

0.00 

0.09 

2.69 

0.384 

0.550 

<0.001* 

0.321* 

WOL 

4.00 

3.65 

0.00 

0.05 

2.08 

0.362 

0.466 

<0.001* 

0.243* 

WYS 

6.00 

5.02 

0.24 

0.43 

3.00 

0.488 

0.599 

<0.001 

0.204* 

SUW 

5.80 

4.86 

0.20 

0.38 

2.93 

0.600 

0.608 

0.083 

0.034 

ZAW 

4.60 

4.09 

0.00 

0.10 

2.42 

0.392 

0.512 

<0.001* 

0.254* 

ITA 

3.80 

3.41 

0.20 

0.14 

2.06 

0.430 

0.471 

0.172 

0.111 

JAP 

4.20 

4.20 

0.73 

1.25 

2.46 

0.414 

0.546 

0.050 

0.276 


collected the specimens in wastelands or meadows, as is typical 
for the species. 

Sample collection. -We collected 184 specimens of A. 
bruennichi in July and August. In Poland, approximately 25 
specimens were collected from each location, resulting in 150 
specimens collected in 2012. Additionally, in 2013, we reeeived 
specimens from Italy (n — 20) and Japan (n— 14). The collected 
material was stored in 95% ethanol at the Museum and 
Institute of Zoology (Polish Academy of Scienees), in Warsaw, 
Poland. 

DNA extraction and amplification.—DNA extractions using 
a single leg from each spider were performed using the 
NucleoSpin® Tissue kit (Macherey-Nagel) according to the 
manufacturer’s protocol. We aimed to amplify the following 
nine microsatellite loci in two multiplex reactions, as described 
by Krehenwinkel & Tautz (2013): MA3, MA5, MA7, MAI 1, 
MAI3, MA27, MA33, MA35 and MA57. Each forward 
primer was labelled with one of the following fluorescent dyes 
(WellRED Dyes): Dye2, Dye3 or Dye4. The reaction mixture 
contained: 1.5 pi of the mixture of primers (“forward” and 
“reverse” for each locus, each primer at 2 pmol/pl); 7.5 pi PCR 
Master Mix (QIAGEN); and 1 pi of DNA extract. The 
reaction mix was then filled to 15 pi with water for PCR 
(Sigma-Aldrich). The reactions were performed under the 
following conditions: 15 min at 95°C; followed by 40 cycles of: 
30 s at 94°C, 90 s at 57°C, 90 s at 72°C; with a final 30 s at 
94°C, 90 s at 57°C, and 10 min at 72°C. The genotyping 
analyses were performed using a CEQ™ 8000 sequencer 
(Beckman Coulter). 

Table 3.—Genetic differentiation among populations of A. 
bruennichi , estimated as F S t- (* — statistically significant values after 
Bonferroni correction, 560 randomizations, adjusted P-valuc at y. = 
0.05 was 0.002). The overall F s T = 0.058 (95% Cl 0.023 0.118). 


Population 

TOR 

WOL 

WYS 

SUW 

ZAW 

ITA 

JAP 

WAR 

0.010 

0.005 

-0.003 

0.022 

0.007 

0.041 

0.182* 

TOR 


0.024 

-0.010 

0.020 

0.008 

0.063* 

0.239* 

WOL 



0.024 

0.049 

0.002 

0.033 

0.252* 

WYS 




0.002 

0.004 

0.042* 

0.206* 

SUW 





0.030 

0.031* 

0.163* 

ZAW 






0.023* 

0.243* 

ITA 







0.238* 


Data analyses. —For each locus within each population, 
deviation from Hardy-Weinberg equilibrium (hereafter HWE) 
and linkage disequilibrium (LZ>) was assessed using Fisher’s 
exact test in Genepop version 4 (Raymond & Rousset 1995; 
Rousset 2008), with the following settings: 10,000 dememori¬ 
zation steps, followed by 1,000 batches and 10,000 iterations. 
Basic genetic indices were calculated for each locus within each 
population as follows: (i) the number of alleles (A); (ii) the 
allelic riehness, or the number of alleles corrected for sample 
size using the rarefaction method with a sample of 14 
individuals ( R A ; Petit et al. 1998); (iii) the effeetive number 
of alleles ( E A ); (iv) the number of private alleles ( R P ); (v) the 
observed (H 0 ) and expected heterozygosity (ff E ) (Nei & 
Roychoudhury 1974); and (vi) the inbreeding coefficient (F 1S ). 
These analyses were performed using the programs GenAlEx 
(Paekall & Smouse 2001), FSTAT version 2.9.3.2 (Goudet 
2001), and HP-RARE (Kalinowski 2005). Moreover, for each 
population we estimated mean values of A, R A , Rp, Hq H e 
and Fi s . 

Genetic differentiation among populations was assessed as 
Fsr (Weir & Coekerham 1984). Pairwise Fsj values and their 
significance, as well as overall F s t with corresponding 95% 
confidence intervals were calculated in FSTAT. 

The Bayesian-elustering method STRUCTURE version 
2.3.4 (Pritchard et al. 2000) was used to examine how well 
the predefined “populations” corresponded to genetic groups 
(K). STRUCTURE was run 10 times for each user-defined K, 
with an initial burn-in of 500,000, and 1,000,000 iterations of 
the total data set. The admixture model of ancestry and the 
correlated model of allele frequencies were used. Sampling 
location was not used as prior information. Next, we 
examined AK statistics to identify the largest change in the 
estimates of K produced by STRUCTURE, as AK may 
provide a more realistic estimation of K than those based on 
likelihood (Evanno et al. 2005). STRUCTURE was run for 
user defined values of K = 1-8. 

To obtain the AK, we used STRUCTURE HARVESTER 
version 0.6.94 (Earl & von Holdt 2011). We then applied 
CLUMPP version 1.1.2 (Jakobsson & Rosenberg 2007) to 
average the multiple runs given by STRUCTURE and correct 
for label switching. The output from CLUMPP was visualized 
with DISTRUCT version LI (Rosenberg 2004). 

As the identification of genetic structure in STRUCTURE 
relies on HWE optimization, and the majority of the sampling 
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Figure 2.—Genetic structure in eight studied populations of A. bruennkhi , inferred using STRUCTURE: A. bar plots for four genetic clusters; 
B. estimated mean (and SD) likelihoods; C. A K curves as a function of K. In bar plots each individual is represented by a vertical bar partitioned 
into segments. The length of each segment describes the estimated membership proportions to each of the genetic clusters. A K suggested the 
presence of four genetic clusters, however the mean likelihood was the highest and had the smallest variance for K = 1. 


sites in our study were not in HWE , we also obtained an 
additional representation of the genetic structure using 
principal component analysis (PCA). This multivariate de¬ 
scriptive method is not dependent on any model assumption 
and ean thus provide a useful validation of the Bayesian 
clustering output (Patterson et al. 2006; McVean 2009; 
Francois & Durand 2010). We used the R package ADEGE- 
NET version L3.4 (Jombart 2008) to carry out standard PCA. 
The results of the analysis were presented graphically along 
the first and second axes according to the highest Eigen values. 

RESULTS 

We failed to amplify the loci MA11 and MAI3 in the 
majority of samples from Poland. Locus MA27 was mono- 
morphic, while MA57 failed to amplify in samples from Italy. 
Hence, final analyses were performed for five microsatellite 
loci, i.e., MA3, MAS, MA7, MA33 and MA35. 

The most polymorphic locus at the population level was 
MA3 (A ranging from 6 to 11 alleles, R from 5.4-8.7), except 
for the JAP population, where M A35 had the highest number 
of alleles (A = 6). In the remaining loci, we observed rather low 
levels of polymorphism, usually from 2 to 5 alleles (Table 1). 
There was no significant linkage disequilibrium (LD) among 
loci within populations (1,600 permutations, adjusted P-value 
after Bonferroni correction: 0.006). In the majority of 
populations, two out of five loci were not in HWE (Table 


1) . Exceptions were populations SUW and ITA, for which 
only one locus exhibited a deviation from HWE. In all cases, 
significant deviation from HWE was due to heterozygote 
deficiency. The inbreeding coefficient (Fi S ) was significant for 
locus MA33 in five of the eight investigated populations. 
Other than MA33, a significant F 1S was also found at locus 
MA3 in TOR, and MA7 in ITA (Table 1). In general, the 
observed heterozygosity was greatest at loci MA3 and MAS in 
populations from Poland and ITA, while in JAP the highest 
H 0 was found at locus MA35 (Table 1), 

For mean indicators of genetic diversity (Table 2), 
populations could be divided into a first group with A around 
6.0 and R > 4.60 (WAR, WYS, SUW), and a second group 
consisting of all remaining populations (A < 4.60; R < 4.09). 
In general, allelie diversity and allelic richness were greater in 
populations from Poland, as opposed to Italy or Japan (Table 

2) . A similar pattern was observed for H 0 . However, the 
highest mean number of private alleles (Pa) and the greatest 
value for private allelic richness (P P ) were observed in the case 
of JAP. 

Five out of eight investigated populations were not in HWE ; 
in four of these F IS was found to be signifieantly higher than 
zero (Table 2). Only one population from Poland (SUW), 
seemed to be in HWE , along with ITA and JAP. 

Genetic differentiation (Fst) was very limited and non¬ 
significant among the populations from Poland (Table 3). We 
found small but significant F§ T values in pairwise comparisons 
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Figure 3.—Principal Component Analysis of A. bruennichi genotypes in eight examined populations. The results are presented along the first 
and seeond axes according to the highest eigenvalues. 1 WAR; 2 - TOR; 3 WOL; 4 WYS; 5 SUW; 6 ZAW; 7 ITA; 8 JAP. 


between populations from Poland and Italy. Genetic differ¬ 
entiation between JAP and the remaining populations was 
high and significant (Table 3). 

Results from STRUCTURE analyses indicated the highest 
mean likelihood for K= 1 and low variance among iterations, 
however a high likelihood was also found for K = 4 (Fig. 2A, 
B). Individuals from Poland and Italy were not assigned to a 
particular cluster, but rather had a mixed probability of 
ancestry from four genetic groups. The population from Japan 
consisted of individuals from a single cluster. The A K 
supported division of the observed genetic variability into 
four genetic clusters (Fig. 2C), and PCA analysis in turn 
suggested a division of populations into two genetic groups 
(Fig. 3), of which the first consisted of the populations from 
Poland and Italy and the second (mainly) of samples from 
Japan. 

DISCUSSION 

Our study indicated that genetic differentiation among 
populations of A. bruennichi in Poland was very low. Counter 
to our expectations, we infer that the strong dispersal abilities 


of this species ensure intensive gene flow, resulting in 
homogeneity of all the Polish populations, which are probably 
still undergoing range expansion. Similar results have been 
obtained for other populations of orb-weaving spiders 
(Ramirez & Fandino 1996; Ramirez & Haakonsen 1999; Lee 
et al. 2004; Jung et al. 2006), and indeed strong gene flow is 
frequently observed between populations of species with high 
dispersal abilities (Waples 1987; Bohonak 1999). However, 
significant genetic differences for species with high dispersal 
abilities have also been recorded (Sielezniew & Rutkowski 
2012; Lemic et al. 2013; Pentek-Zakar 2015). 

According to Vucetich & Waite (2003), populations at the 
edge of a species’ range are generally smaller than those at the 
core, and observations for A. bruennichi suggest that the 
north-eastern (SUW) population is indeed smaller than others 
(e.g., 0.2 ind. /m 2 in SUW and 1.75 ind. /m 2 in HOR) (W. 
Wawer, unpubl. data). However, neither the new (SUW) nor 
the old (e.g., WYS and WOL) populations are in HWE (Table 
1). In fact, due to a heterozygote deficiency, many populations 
of expansive species are not found to indicate HWE , e.g., 
among Lepidoptera (Wei et al. 2013), Coleoptera (Wu et al. 
2016) and Diptera (Yan et al. 2015). In the case of spiders, 
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HWE has usually been noted in stable-range populations 
(Ramirez & Chi 2004; Ramirez et al. 2007, 2013). 

Allelic richness was found to be higher in Poland than in 
Italy. This observation supported previous results, indicating 
high genetic diversity in expanding populations of the speeies 
(Krehenwinkel & Tautz 2013). It is possible that there is still 
ongoing gene flow among populations in Poland, eliminating 
founder effects and isolation in newly colonized areas 
(Greenbaum et al. 2014). Higher genetic diversity in Poland 
as compared to Italy has also been recorded in other 
invertebrates (Sielezniew et al. 2015), and was explained by 
“the rear edge hypothesis”, whereby rear edge populations are 
usually small in size and are thus characterized by low genetie 
diversity and high inter-population genetic differentiation 
(Hampe & Petit 2005). Furthermore, marked genetic differ¬ 
entiation between Asian and European populations was also 
revealed in our study. Similar results were obtained by 
Krehenwinkel & Tautz (2013), who confirmed genetic 
distinctness when East Asian and Western Palearctic popula¬ 
tions were compared. As the peripheral and recently settled 
(SUV/) population was characterized by the highest value for 
heterozygosity ( H 0 = 0.600) and the lowest inbreeding 
coefficient (Fis = 0.034), it is possible to conclude that a small 
population of A, bruennichi at the edge of its range is 
connected with others by gene flow. Moreover, despite their 
relatively long-term persistence (of ca. 80 years), the remaining 
Polish populations appear to still be in the expansion phase, 
given the lack of HWE. Similar results were obtained by 
Wellenreuther et al. (2011) for populations of the highly vagile 
blue-tailed damselfly, Ischnura eiegans (Vander Linden, 1820). 
Generally, populations displaying a high expansion rate are 
associated by gene flow, such that even long distances or 
geographical barriers do not inhibit the process (Well¬ 
enreuther et al. 2011) And although genetic diversity is 
greater far from the center, in most cases distant populations 
are seen to be similar genetically (Eckert et al. 2008). 

At this stage of the research, we cannot confirm the 
hypothesis of a colonization of Polish territory from both 
the east and west (Urbahski 1948; Barabasz & Gorz 1998). 
Analyses of genetic structure indicated similarity between the 
Polish and Italian populations, and while the Alps are known 
to have aeted as an initial barrier to the possible northward 
expansion of Mediterranean genetic pools (Bilton et al. 1998; 
Hewitt 1999), the Appenine Peninsula might be a possible 
source area of invertebrate populations in Central and 
Northern Europe (Patricelli et al. 2013). This observation 
supports previous results obtained by Krehenwinkel & Tautz 
(2013), suggesting that northern territories of Europe were 
colonized from southern localities rather than from eastern 
areas. 
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